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The  potential  and  current  distributions  are 
calculated  across  the  width  of  a  quantum 
Hall  effect  sample  for  applied  currents  be¬ 
tween  0  p-A  and  225  |jlA.  For  the  first 
time,  both  a  confining  potential  and  a  cur- 
rent-induced  charge-redistribution  potential 
are  used.  The  confining  potential  has  a 
parabolic  shape,  and  the  charge-redistribu¬ 
tion  potential  is  logarithmic.  The  solution 
for  the  sum  of  the  two  types  of  potentials 
is  unique  at  each  current,  with  no  free 
parameters.  For  example,  the  charge-deple¬ 
tion  width  of  the  confining  potential  is  de¬ 
termined  from  a  localization  experiment  by 
Choi,  Tsui,  and  Alavi,  and  the  spatial  extent 
of  the  conducting  two-dimensional  electron 
gas  across  the  sample  width  is  obtained 
from  the  maximum  electric  field  deduced 
from  a  high-current  breakdown  experiment 
by  Cage  and  Lavine,  and  from  the  quantum 


Hall  voltage.  The  spatial  extent  has  realistic 
cut-off  values  at  the  sample  sides;  e.g.,  no 
current  flows  within  55  magnetic  lengths  of 
the  sides  for  currents  less  than  215  |jlA. 

The  calculated  potential  distributions  are  in 
excellent  agreement  with  contactless  elec¬ 
tro-optic  effect  laser  beam  measurements  of 
Fontein  et  al. 
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1.  Introduction 

The  potential  and  current  distributions  within  quantum 
Hall  samples  are  important  aspects  of  the  integer  quan¬ 
tum  Hall  effect  [1-3].  In  this  effect  the  Hall  resistance, 
Rh,  of  the  ith  plateau  of  a  fully  quantized  two-dimen¬ 
sional  electron  gas  (2DEG)  has  the  value  Ri{(i)=h/(e^i), 
where  h  is  the  Planck  constant,  e  is  the  elementary 
charge,  and  i  is  an  integer.  Early  attempts  to  measure 
potential  distributions  across  samples  [4-7]  used  electri¬ 
cal  contacts  to  the  two-dimensional  gas  that  were  placed 
within  the  sample  interior.  The  potentials  were  found  to 
vary  throughout  the  entire  sample.  There  was  concern, 
however,  that  the  electrical  contacts  themselves  signifi¬ 
cantly  altered  the  potential  distributions.  Eontein  et  al. 
[8]  have  made  contactless  measurements  of  potential 


distributions  using  a  laser  beam  and  the  electro-optic 
Pockels  effect.  They  observed  major  fractions  of  the 
quantum  Hall  voltage  occurring  near  the  sides  of  the 
sample,  but  also  significant  contributions  within  the 
interior.  Valid  predictions  of  the  potential  distribution 
across  quantum  Hall  samples  should  agree  with  their 
results. 

In  this  paper  we  calculate  the  potential  distributions 
across  the  sample  for  applied  currents  /sd  between  0  pA 
and  225  pA  by:  (a)  assuming  a  parabolic  confining 
potential  for  the  charge  carriers  and  using  parameters  of 
the  parabola  obtained  experimentally  by  Choi,  Tsui,  and 
Alavi  [9];  (b)  assuming  an  applied  current-induced  loga¬ 
rithmic  charge-redistribution  potential  for  the  charge 
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carriers  that  is  similar  to  that  of  Balaban,  Meirav,  and 
Shtrikman  [10],  but  with  very  different  cut-off  values 
for  the  spatial  extent  of  the  potential;  (c)  assuming  that 
the  width  of  the  conducting  region  varies  with  applied 
current  because  a  Lorentz  force  deflects  the  conducting 
electrons  slightly  towards  one  side  of  the  sample;  (d) 
using  the  maximum  electric  field  deduced  by  Cage  and 
Lavine  [11]  from  a  breakdown  experiment  at  high  cur¬ 
rents  to  determine  the  cut-off  value  for  the  spatial  extent 
on  one  side  of  the  sample;  and  (e)  using  the  quantum 
Hall  voltage  to  determine  the  cut-off  value  of  the  spatial 
extent  on  the  other  side  of  the  sample.  The  calculated 
potential  distributions  are  in  excellent  agreement  with 
the  measurements  of  Fontein  et  al.  [8]. 

2.  Coordinate  System 

The  coordinate  system  is  shown  in  Fig.  1.  For  conve¬ 
nience  in  writing  the  equations,  the  origin  is  located  at 
the  source  S  and  is  halfway  across  the  sample  width  w. 
The  sample  labeling  is  chosen  to  be  consistent  with 
previous  work  [11].  Potential  probes  1  and  2,  and  the 
drain  D,  are  not  shown.  The  positive  x  axis  points  along 
the  sample  in  the  general  direction  of  the  externally 
applied  current  /so.  The  positive  y  axis  is  chosen  as 
indicated.  Therefore  the  positive  z  axis  points  downward 
for  a  right-handed  coordinate  system.  The  magnetic 
field  B  also  points  downward,  simply  to  be  consistent 
with  results  from  the  breakdown  experiment  [11]  that 
will  be  used  in  Secs.  4  and  5. 

Note  that  the  conducting  charges  are  electrons  with 
charge  q  =  -  e .  This  is  taken  into  account  throughout  the 
paper;  it  is  necessary  to  do  so  because  the  signs  of  both 
the  confining  potential  and  the  charge-redistribution 
potential  depend  on  the  sign  of  the  charge  carriers. 
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Fig.  1.  The  coordinate  system  used  in  this  paper.  Only  part  of  the 
sample  is  shown.  The  origin  is  located  at  the  source  S,  halfway  across 
the  sample  width  w.  The  dotted  curves  indicate  the  electron  flow 
pattern  for  this  magnetic  field  direction.  Fl  is  the  Lorentz  force  on  the 
conducting  electrons  and  Fc  is  the  Coulomb  force.  B  is  the  magnetic 
field,  V  is  the  electron  velocity,  and  Va  and  Vb  are  the  potentials  on 
either  side  of  the  sample. 


The  Lorentz  force  Fi^  =  ev  X  B  is  in  the  positive  y 
direction.  This  force  deflects  the  electrons  slightly  to  the 
right  until  it  is  matched  by  the  Coulomb  repulsive  force 
Fc  =  -eE  [12].  A  charge-redistribution  of  the  2DEG 
results  from  this  deflection.  Also  because  of  the  Lorentz 
and  Coulomb  forces,  the  electrons  enter  the  corner  of  the 
source  at  y  =  -  w/2  for  this  magnetic  field  direction  and 
exit  at  the  corner  h-w/2  of  the  drain — in  agreement  with 
the  experiment  of  Klass  et  al.  [13].  We  assume  that  the 
electrons  spread  out  across  the  sample  interior  in  agree¬ 
ment  with  the  experiment  of  Fontein  et  al.  [8].  Potential 
probes  4  and  6  are  near  the  potential  of  the  source. 
Probes  3  and  5  are  near  the  potential  of  the  drain,  and 
have  a  positive  potential  relative  to  the  source  for  these 
current  and  magnetic  field  directions.  The  chemical  po¬ 
tential  (pA  =  Va  is  therefore  positive  relative  to  the  chem¬ 
ical  potential  (^  =  Vb  on  the  opposite  side  of  the  sample. 

3.  Confining  Potential 

We  begin  the  calculations  with  a  confining  potential 
to  prevent  the  2DEG  from  spilling  out  the  sides  of  sam¬ 
ples.  Choi,  Tsui,  and  Alavi  [9]  performed  an  experiment 
on  mesa-etched  GaAs/AlGaAs  heterostructure  samples 
in  zero  magnetic  field.  They  then  used  one-dimensional 
localization  theory  to  deduce  the  charge-depletion 
widths.  A,  of  the  confining  potentials,  and  found  that  A 
was  (0.5  ±  0.2)  p,m  for  a  2DEG  of  surface  number 
density  Ws  =  i{eBlh)  =  1.5  X  lO'Vcm^.  We  will  use  their 
results  to  define  the  depletion  width  of  the  confining 
potential  for  a  mesa-etched  sample. 

3.1  Charge-Depletion  Region 

Eigure  2  (a)  shows  a  schematic  of  the  charge  distribu¬ 
tion  in  the  GaAs/AlGaAs  interface  region  near  one  side 
of  the  mesa  when  there  is  no  applied  magnetic  field.  The 
GaAs  layer  of  our  sample  [11]  has  a  residual  donor 
density  of  about  1  X  lO'Vcm^,  while  the  donor  concen¬ 
tration  in  the  AlGaAs  layer  is  about  1  X  lO'Vcm^  and 
Hs  =  5.94  X  lO'Vcm^.  There  is  an  ionized  donor  atom  in 
the  AlGaAs  layer  for  every  electron  in  the  2DEG  but, 
unlike  Choi,  Tsui,  and  Alavi  [9],  we  assume  the  ionized 
donor  atoms  are  distributed  over  a  volume  rather  than  in 
a  surface  sheet  with  density  n^.  The  confining  potential 
is  generated  from  electron  surface  charges  on  the  side  of 
the  mesa,  as  indicated  in  the  figure.  There  is  an  ionized 
donor  atom  or  ionized  impurity  site  in  the  charge-deple¬ 
tion  region  for  every  surface  charge. 

We  assume  a  homogeneous  charge-depletion  region 
in  Eig.  2  (b).  The  depletion  width  A  for  a  homogeneous 
three-dimensional  material  is  [14] 

A  =  {2eA^JeNuy'\  (1) 
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Fig.  2.  (a)  Schematic  diagram  of  the  GaAs/ AIGaAs  interface  region 

near  one  side  of  the  mesa.  See  Sec.  3.1  for  further  explanation,  (b)  The 
ionized  donor  charge-depletion  density  distribution  Ad-  (c)  The  con¬ 
fining  potential  Vc  for  negatively  charged  surface  states. 


where  Ss  =  kso  is  the  dielectric  permittivity  of  the  semi¬ 
conductor,  K  is  the  dielectric  constant  (k  =  13.1  for 
GaAs),  So  is  the  permittivity  of  vacuum,  Vm  is  the  value 
of  the  confining  potential  at  ±w/2  and  Ad  is  the  average 
density  of  ionized  donors  and  impurity  sites  in  the 
charge-depletion  region.  We  selected  the  value  of  the 
charge-depletion  width  to  be  Zl  =  0.5  p.m  [15].  This 
value  is  consistent  with  the  results  of  Choi,  Tsui,  and 
Alavi  [9].  We  chose  the  value  of  Vm  to  be  one-half  the 
1.50  V  separation  between  the  valence  and  conduction 
bands  of  GaAs  at  1  K  [14],  or  =  0.75  V.  The  value  of 
the  average  charge-depletion  density  from  Eq.  (1)  is  thus 
Ad  =  4.3  X  lO'Vcm^  which  seems  quite  reasonable. 

3.2  Confining  Potential  Equation 

A  homogeneous  charge-depletion  region  results  in  a 
parabolic  confining  potential  y,,  with  the  origin  at 
y  =  A  =  w/2-A,  as  indicated  schematically  in  Fig.  2  (c). 
The  confining  potential  is  negative  because  the  charges 
on  the  side  of  the  mesa  are  electrons. 

The  equations  for  the  confining  potential  Vc  and  its 
electric  field  Ec  =  -  VVc  are 


Tc(y)  =  -  a(y+Xf  and  Edy)  =  2a(y+X)  (2c) 

for-y  <-y<-A, 


where  a  =  VJA^  =  3.0  X  10'^  Vlwd  for  A  =  0.5  pm  and 
Vm  =  0.75  V,  and 


A 


A  . 


(3) 


3.3  Confining  Potential  at  Isd  =  0  /tA 

Given  the  values  of  A  and  Vm,  there  is  a  surprising 
amount  that  can  be  deduced  about  the  electron  states  of 
the  confining  potential  when  the  magnetic  field  is  ad¬ 
justed  to  be  at  the  center  of  the  i  =  2  quantum  Hall 
plateau  and  Ad  =  0  pA.  Since  there  is  no  applied  current, 
and  therefore  no  Hall  voltage,  the  Fermi  energy  sp  is 
constant  across  the  sample  width  and  is  located  halfway 
between  Fandau  levels.  Under  these  conditions,  states  of 
the  lowest  (A  =  0)  Fandau  level  are  occupied  up  to  the 
Fermi  energy  sp  =  ^u)J2,  no  states  are  occupied  in  the 
second  (A=  1)  Fandau  level,  a>c  =  eB/m*  is  the 
cyclotron  angular  frequency,  m*  is  the  reduced  mass  of 
the  electron  (0.068  times  the  free  electron  mass  in 
GaAs),  and  ^  =  hl2'TT.  References  [11,16-18]  describe 
how  these  states  can  be  defined  in  the  Fandau  gauge. 

Figure  3  shows  a  schematic  drawing  of  the  energy  of 
the  confining  potential  for  Ad  =  0  pA,  with  greatly 
exaggerated  values  of  A  and  fiwc,  and  only  a  small  frac¬ 
tion  of  allowed  states.  The  occupied/unoccupied  states 
are  indicated  as  solid/open  circles,  and  the  occupied 
(filled)  states  lie  between  y^ax  and  ymin  =  -ymax.  In  the 
presence  of  the  magnetic  field,  electrons  of  the  2DEG 
occupy  Fandau  level  states  that  penetrate  into  the 
charge-depletion  regions  near  the  mesa  edge,  and  cur¬ 
rent  circulates  around  the  sample  periphery.  Under  these 
conditions 

/io) 

^c(ymax)  —  2  —  ^^cCypiax)  —  (ymax 

=  (ymax- A)^  ,  (4) 
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where  A  =  w/2  -  A. 

The  occupied  states  of  the  right-hand  side  (rhs)  con¬ 
fining  potential  generate  a  total  current  Ic  (rhs)  that  is 


Tmax 

/c(rhs)  =  j  Jc(y)dy 
A 


Tmax 

j  (T^Ec(y)dy  = 

A 


-  ^  -  Vc(A)],  =  -  ,  (5) 

fin  fiH 


where  Jc(y)  is  the  current  density,  a-j^y  is  the  off-diagonal 
conductivity  tensor  component,  VcCymax)  =  -a(ymax- 
A)^,  and  VcCA)  =  0.  In  the  absence  of  significant  dissipa¬ 
tive  scattering  on  the  quantum  Hall  plateau,  cr^y  =  1/Rh 
[12].  Similarly, 


-A 


/c  (Ihs)  = 


a^EcCy)  dy  = 


Tmin 


y  min 


Fig.  3.  Schematic  drawing  of  the  energy  of  the  confining  potential 
Vc  across  the  sample  when  /so  =  0  fxA.  Values  of  the  charge-depletion 
width  A  and  the  Landau  energy  level  spacing  hiOc  are  greatly  exagger¬ 
ated.  The  occupied/unoccupied  states  of  the  first  two  Landau  levels  are 
shown  as  solid/open  circles.  The  occupied  (filled)  states  lie  between 
the  locations  y^ax  =  -  ymin- 


4.  Charge-Redistribution  Potential 


-  ^  [Vc(-A)  -  yc(y„,i„)],  =  -  ,  (6) 


where  VcCymm)  =  -  fl(ymm+A)l 

It  follows  from  Eqs.  (2)  to  (6)  for  the  12  906.4  U, 
i  =  2  plateau  at  12.3  T,  for  the  400  p.m  wide  sample  of 
Ref.  [11],  and  for  /sd  =  0  p.A  that 

'■  fr"**  =  It; 


Section  2  noted  that  the  Lorentz  force  exerted  on  the 
conducting  electrons  causes  deviations  -e8cr(y)  from 
the  average  surface  charge  density  -  ecr^e  =  -en^  = 
-  ie^B/h  of  the  2DEG  charge-redistribution  across  the 
sample  width.  The  resulting  charge-redistribution 
potential,  Vr(y),  arising  from  applied  currents  would  be 
a  linear  function  of  y  if  the  mobile  electrons  occupied  a 
three-dimensional  volume.  They  occupy  a  two-dimen¬ 
sional  sheet,  however,  and  MacDonald,  Rice,  and 
Brinkman  [19]  expressed  this  charge-redistribution  self- 
consistently  in  terms  of  a  charge-redistribution  potential 
as 


and 


ymax  =  -  ymin  =  199.559  p,m  ,  (8) 


Vr(y)  =  - 


ItTKBo 


8(T(y')  In 


y  -y 


dy'  ,  (10) 


w 

2 


^max 


=  0.441  p,m. 


(9) 


Thus,  a  rather  large  0.8 1  p-A  current  circulates  around 
the  sample  at  12.3  T  when  /sd  =  0  pA,  A  =  0.5  pm  and 
Vm  =  0.75  V.  The  maximum  extent  of  this  current  is  60 
times  farther  from  the  sides  of  the  sample  than  that 
produced  by  skipping  orbits  bouncing  off  of  a  hard  wall 
with  a  cyclotron  radius  or  magnet  length  Zg  =  (h/eBY'^ 
of  7.3  nm. 


where 


8o-(y)  = 


im  *  d^ 
~hB  d/ 


^  d/ 


yr(y), 


(11) 


as  shown  in  Appendix  A.  Riess  [20]  extended  this  po¬ 
tential  to  a  2DEG  with  finite  thickness.  Thouless  [21] 
then  found  an  analytic  logarithmic  approximation  of  this 
potential  far  from  the  sample  sides,  and  Beenakker  and 
van  Houten  [22]  then  approximated  the  near-edge  be¬ 
havior  by  introducing  a  cut-off  at  a  distance  ^  from  the 


532 


Volume  100,  Number  5,  September-October  1995 

Journal  of  Research  of  the  National  Institute  of  Standards  and  Technology 


sample  side,  and  a  linear  extrapolation  for 


Vr(y)  =  -^  I  ^VAy')ln 


-ly'-yl 

w 


dy'  (12) 


from  ly  I  =  w/2  -  ^  to  ly  I  =  w/2.  The  characteristic  length 
^is  illlrra*  =  ieVilTTKSohw^)  ^  where  Zg  = 
is  the  magnetic  length  and  a*  =  4TTKSo^^lm*e^  is  the 
effective  Bohr  radius  in  SI  units.  Our  values  of  Zb  , 
and  fl*  for  the  i  =  2  plateau  at  12.3  T  are  3.3  nm,  7.3  nm, 
and  10.2  nm,  respectively. 

Balaban,  Meirav,  and  Shtrikman  [10]  used  a  nonlin¬ 
ear  (quadratic)  extrapolation  near  the  sample  sides  and 
obtained  the  charge-redistribution  potential 


the  /sD  =  0  pA  conditions  of  Sec.  3.3  and  then  reach 
the  electric  field  values  necessary  for  quasi-elastic  inter- 
Landau  level  scattering  (QUILLS)  transitions  [11,  16- 
18,25-28]  at  high  currents. 

We  use  the  same  form  for  the  charge-redistribution 
potential  as  Balaban  et  al.  [10],  but  with  a  different 
geometrical  factor  and  very  different  cut-off  values,  y^in 
and  ymax,  which  vary  with  applied  current.  Our  charge- 
redistribution  potential  is 


44h 

ymax  H-  w/2 

-1 

In 

z. 

+  w/2 

2 

_  w/2  -  ymax_ 

Z 

-  w/2 

(15) 


for 


w  w 

-  2  <  I'm"  -  y  -  Jmax  <  2 


Vr(y)  =  - 


7si>4h 

w  ^  +  5 
In  ^  H- 

-1 

In 

y  +  w/2 

2 

8  2^  _ 

y  -  w/2 

(13) 


where 


4  =  /sD  -  4  (rhs)  -  4  (Ihs). 


(16) 


for  ly  I  <  w/2  -  S,  where  8  =  Ib  for  the  i  =  2  plateau,  and 
8  is  not  the  differential  8  of  Eq.  (10).  They  successfully 
used  this  potential  to  describe  the  sample-width  depen¬ 
dence  for  breakdown  at  small  currents,  but  could  not 
account  for  the  larger  breakdown  currents  observed  in 
other  experiments  [11,  23-28].  Their  geometry  factor  is 


g+  g]-' 
2^  . 


0.08 


(14) 


for  our  values  of  ^  and  S  at  w  =  400  pm. 


4.1  Charge-Redistribution  Potential  Equation 

The  charge-redistribution  potential  described  by 
Eq.  (13)  was  calculated  for  an  infinite  square-well  con¬ 
fining  potential,  and  must  be  modified  for  use  with  a 
more  realistic  confining  potential.  To  do  this  correctly 
would  require  a  numerical  solution  of  Eq.  (12),  with  the 
confining  potential  included,  as  is  discussed  in 
Appendix  A.  We  approximated  this  numerical  solution 
(and  then  tested  the  approximation)  by  using  the  form  of 
the  potential  in  Eq.  (13)  but  introducing  two  parameters, 
ymin  and  y^ax,  that  alter  the  charge-redistribution  poten¬ 
tial  due  to  the  presence  of  the  quadratic  confining 
potential. 

It  was  necessary  to  do  this  because  the  potential 
distribution  of  Eq.  (13),  with  a  cut-off  distance  S  =  Zg, 
gave  the  correct  quantum  Hall  voltage  Vh  =  4h/sd  across 
the  sample,  but  the  electric  field  Ey  =-  VW  did  not 
increase  quickly  enough  for  increasing  current  to  satisfy 


4(rhs)  and  4(lhs)  are  defined  by  Eqs.  (5)  and  (6),  and 
the  geometry  factor  G  in  Eq.  (15)  is 


G{W,  ymax)  = 


ymax  +Wl2  ‘ 

w/2  -  ymax  _ 


.  (17) 


We  assume  G  is  current-independent,  and  assign  the 
value 


G  =  0.147  (18) 

to  Eq.  (17)  by  using  the  value  of  ymax  =  199.559  pm 
found  in  Sec.  3.3  for  /sd  =  0  pA  and  w  =  400  pm. 
Our  value  of  G  is  thus  somewhat  larger  than  the  value 
G  =  0.08  that  would  be  used  by  Balaban  et  al.  [10].  The 
cut-off  values 


=  w/2  -  ymax  and  Smln  =  w/2  H-  ymm  (19) 

will  be  determined  in  Sec.  5.  Appendix  B  discusses  the 
agreement  between  our  Eq.  (15)  and  the  self-consistent 
Eqs.  (10)  and  (11). 

The  electric  field  Ey  =-Wy  due  to  redistribution  of 
the  2DEG  with  applied  current  is 


EAy)  =  G 


w 


[{wHy-y^] 


■  (20) 


We  now  have  nearly  all  the  information  necessary  to 
determine  the  potential  and  current  distributions. 
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5.  Calculations 

Figure  4  shows  the  confining  potential  -  Vc(y)  and 
the  charge-redistribution  potential  -  V^iy)  across  the 
sample  for  greatly  exaggerated  values  of  A,  S„ax,  and 
Smin,  and  for  an  arbitrary  value  of  4,  where  4  is  defined 
by  Eq.  (16).  Vr  becomes  infinite  at  ±w/2,  but  that  is  of 
no  concern  because  it  is  only  the  occupied  states  which 
contribute  to  the  Hall  voltage,  and  those  states  occur 
only  between  ymax  and  ymin.  The  potentials  are  therefore 
finite  and  well-behaved  in  the  region  of  interest. 


Fig.  4.  The  confining  potential  -  Vc(y)  and  the  charge-redistribution 
potential  -  VrCv)  across  the  sample  for  greatly  exaggerated  values  of 
4  s,„„  =  w/2-)>max,  and  w/2  +  yrrin. 

5.1  Total  Potential 

Of  course  the  electrical  transport  properties  depend 
on  the  total  potential  VtCy),  but  we  can  unambiguously 
separate  VAy)  into  the  confining  and  charge-redistribu¬ 
tion  potential  components 

Vt(y)  =  VAy)  +  VAy)  ■  (21) 

We  have  uniquely  defined  the  potentials  VcCy)  and  VAy) 
in  Eqs.  (2)  and  (3)  of  Sec.  3.2  and  Eqs.  (15)  to  (18)  in 
Sec.  4.1,  plus  Eqs.  (5)  and  (6)  in  Sec.  3.3.  The  current- 
independent  parameters  for  the  confining  potential  and 
the  charge-redistribution  potential  are:  A  =  0.5  p.m, 
Vm  =  0.75  V,  and  G  =  0.147.  Eor  a  given  sample  we  know 
the  applied  current  /sd  and  the  sample  width  w,  but  there 
are  still  two  free  parameters:  y^ax  and  ymin. 

Ordinarily,  it  would  not  be  possible  to  uniquely  deter¬ 
mine  the  values  of  y^ax  and  y^in  since  the  only  other  piece 
of  information  is  that  the  quantum  Hall  voltage  Vh  is 

Vh  =  4h/sD  =  Vt(ymin)  —  Vt(ymax),  (22) 

and  there  is  a  range  of  values  for  y^ax  that  satisfies  this 
equation.  It  is  possible,  however,  to  determine  the  value 


of  ymax  for  a  particular  type  of  experiment,  and  we 
believe  that  the  results  are  representative  of  most  other 
experiments  since  our  calculations  agree  with  the  ex¬ 
perimental  data  of  Eontein  et  al.  [8].  We  first  note  that 
E(y)  =  -  VV(y).  Therefore 

E't(ymax)  “  4c(yniax)  "b  4'r(yniax)  .  (73) 

In  an  experiment  described  in  Ref.  [11]  we  measured 
the  quantized  longitudinal  voltage  drops  along  a  GaAs/ 
AlGaAs  sample  between  potential  probes  4  and  6  of 
Eig.  1  at  high  currents,  and  deduced  the  maximum  elec¬ 
tric  field  Eniax  from  a  quasi-elastic  inter-Landau  level 
scattering  model.  The  results  were 

Emax  =1.1X10®  V/m  @  /sD  =  215  p.A  (24a) 
and 

4„ax  =  4.2X10®  V/m  @  /sd  =  225  pA.  (24b) 

The  value  Emax  =  1 . 1 X 10®  V/m  at  /sd  =  215  pA  was  just 
sufficient  to  excite  the  lowest,  M  =  1,  QUILLS  transi¬ 
tions  [11,25-28].  It  is  clear  from  Eig.  4  that  Emax  will 
occur  at  y„ax,  so 

£'t(ymax)  =  Emax.  (25) 

We  can  therefore  use  Eqs.  (23)  and  (24)  to  determine 
ymax,  and  then  Eq.  (22)  to  obtain  ymin  for  the  sample  of 
Ref.  [11].  Note  that  changing  the  values  of  ymax  and  ymin 
also  alters  the  values  of  4(rhs),  4(lhs),  and  thereby  the 
value  of  4  in  Eqs.  (5),  (6),  and  (16).  Thus  there  are  no 
free  parameters ,  and  one  can  obtain  unique  solutions  to 
the  total  potential  and  other  transport  properties. 

5.2  Results 

Relevant  values  for  the  solution  at  Ed  =  0  pA  are 
shown  in  Table  1.  Most  were  calculated  in  Sec.  3.3;  the 
remainder  were  found  from  Eqs.  (2),  (3),  and  (15)  to 
(22).  Note  that  ymax  and  ymin  are  predicted  to  be  about  60 
magnetic  lengths  from  the  sides  of  the  sample. 

We  calculate  the  values  shown  in  Table  1  at 
/sd  =  215  pA  by  increasing  the  value  of  ymax  until 
Et(ymax)  =  1.1  X  10®  V/m,  adjusting  the  value  of  ymin  to 
obtain  the  correct  Hall  voltage,  and  remembering  that 
changing  the  values  of  ymax  and  ymin  also  changes  the 
values  of  4(rhs),  4(lhs),  and  4.  The  solution  is  unique, 
with  no  free  parameters.  The  same  procedure  is  done  at 
/sd  =  225  pA,  except  that  the  value  of  ymax  is  increased 
until  E,(ymax)  =  4.2  X  10®  V/m.  Note  in  Table  1  that  ymax 
is  still  about  13  magnetic  lengths  away  from  the  side  of 
the  sample  at  /sd  =  225  pA. 
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Table  1.  Values  of  some  quantities  obtained  from  the  calculations  in  Sec.  5  for  Isd  =  0  jjlA,  25p.A,  215  p-A,  and  225  pA.  The  quantities 
common  to  all  the  calculations  are  i  =  2,B=  12.3  T,  k  =  13.1,  w  =  400  pm,  A  =  0.5  pm,  =  0.75  V,  a  =  3.0X  10‘^  V/m^,  A  =  199.500 
pm,  G  =  0.147,  and  4  =  7.3  nm.  See  Secs.  2-5  for  the  definitions  of  these  symbols 


/sD 

/c(rhs) 

(pA) 

/Plhs) 

(pA) 

I, 

(pA) 

^max 

(pm) 

iVmin 

(pm) 

^cCVmax) 

(V) 

FCymax) 

(V) 

FcCVmax) 

(MV/m) 

£'r(ymax) 

(MV/m) 

^tCVmax) 

(MV/m) 

^max/ /fi 

0 

0.81 

-0.81 

0.00 

199.559 

-199.559 

-0.010 

0.000 

0.354 

0.000 

0.354 

60.3 

25 

0.94 

-0.68 

24.74 

199.564 

-199.554 

-0.012 

-0.160 

0.382 

0.054 

0.436 

59.6 

215 

2.30 

-0.05 

212.75 

199.599 

-199.515 

-0.030 

-1.392 

0.596 

0.504 

1.100 

54.8 

225 

37.36 

-0.00 

187.64 

199.901 

-198.044 

-0.482 

-1.477 

2.405 

1.795 

4.200 

13.5 

We  also  calculate  the  relevant  quantities  at  /sd  =  25 
pA,  which  is  a  current  often  used  in  precision  quantized 
Hall  resistance  measurements.  In  this  case,  however,  we 
do  not  know  the  value  of  fit  (y max),  so  we  use  a  linear 
interpolation  of  the  value  of  jmax  between  its  values  for 
fsD  =  0  p,A  and  215  pA.  The  quantities  shown  in  Table  1 
for  fsD  =  25  pA  are  relatively  insensitive  to  this  choice 
for  ymax- 

5.3  Plots 

We  now  plot  the  potentials,  using  Eqs.  (2),  (3),  (15) 
to  (19),  and  (21).  Figure  5  shows  yc(y)  and  V^iy)  for  the 
parameters  used  in  Table  1  at  /sd  =  215  pA,  except  that 
the  plot  is  between  ±0.99999  w/2  (±  199.998  pm) 
rather  than  jmax  and  y^m  in  order  to  show  the  sharpness 
of  the  confining  potential  and  the  extent  of  the  charge- 
redistribution  potential  at  these  extreme  values  of  y. 
Figure  6  shows  Vi{y)  plotted  between  jmax  and  ymin  using 
the  parameters  in  Table  1  at  /sd  =  215  pA  and  225  pA. 
Other  than  moving  farther  to  the  right,  the  total  potential 
does  not  significantly  change  shape  with  increasing 
current. 

Figure  7  shows  Vi{y)  at  /sd  =  25  pA.  The  shape  of  this 
predicted  potential  is  in  excellent  agreement  with  the 
experimental  measurements  shown  in  Fig.  6  of  Fontein 
et  al.  [8].  It  is  this  agreement  which  provides  the  best 
verification  of  our  results.  The  “linear”  part  of  the  po¬ 
tential  distribution  within  the  sample  interior,  attributed 
in  Ref.  [8]  to  heating  effects  which  cause  =  Vx//sd  to 
increase,  is  accounted  for  by  our  charge-redistribution 
potential  in  a  sample  which  has  minimal  heating  at  these 
currents  [24]. 

The  electric  fields  E^{y)  =  -  Vyc(y)  and  £’r(y)  =  - 
Vyr(y)  are  shown  in  Fig.  8  for  /sd  =  215  pA;  they  were 
determined  from  Fqs.  (2),  (3),  (18),  and  (20).  The  value 
of  ymax=  199.599  pm  is  such  that  £'t(ymax)  =  1.1X10*’ 
V/m  in  equation  (23).  The  contribution  to  the  total  elec¬ 
tric  field  at  jmax  is  slightly  more  for  the  confining 
potential  than  for  the  charge-redistribution  potential  at 


this  current.  Table  1  shows  that  the  confining  potential 
also  provides  the  dominant  contribution  to  /St(ymax)  at 
other  currents. 


Fig.  5.  yCr)  and  VAy)  plotted  between  ±0.99999  w/2  for  the 
parameters  used  in  Table  1  at  /sd  =  215  p-A.  The  parameters  common 
to  all  plots  in  Figs.  5-1 1  are  i  =  2  (12  906.4  ft),  B  =  12.3  T,  w  =  400 
pm,  K  =  13.1,  A  =  0.5  pm,  y„  =  0.75  V,  and  G  =  0.147. 


The  location,  jmax,  of  the  last-filled  state  on  the  right- 
hand  side  of  the  sample  increases  with  applied  current 
/sd.  We  can  use  Eq.  (A-3)  and  Table  1  to  determine  what 
part  of  this  increase  in  y^ax  is  due  to  the  increase  in  the 
total  electric  field  at  ymax-  The  percentage  contributions, 
relative  to  the  values  of  y^ax  and  £'t(ymax)  at  /sd  =  0  pA, 
are  4  %,  5  %,  and  3  %,  for  /sd  =  25  pA,  215  pA,  and 
225  pA,  respectively.  Therefore,  most  of  the  increase  in 
ymax  is  due  to  the  Forentz  force  pushing  the  electrons 
closer  to  the  side  of  the  sample. 

The  current  density  Jt(y)  for  electrons  moving  in  the 
positive  X  direction  is 

My)  =  EM)  =  ^  [£c(y )  ±  EM)l  (26) 
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Fig.  6.  Vt(y)  plotted  between  ymax  and  >'niin,  using  the  parameters  in 
Table  1  for  /so  =  215  |jlA  and  225  p-A.  The  values  of  ^max  and  jmin  are 
199.599  |JLm  and -199.5 15  [xm,  and  199.901  ixm  and -198.044  jxmfor 
IsD  =  215  fxA  and  225  |xA,  respectively. 


y  (lim) 


Fig.  7.  VtCy)  at  /sD  =  25  |jlA.  This  potential  is  in  excellent  agreement 
with  the  experimental  measurements  shown  in  Fig.  6  of  Fontein  et  al. 
[8]. 


Figure  9  shows  Jt(y)  for  /so  =  25  |xA,  215  ixA,  and  225 
piA.  The  maximum  two-dimensional  current  density  is 
at  jmax,  and  is  85  A/m  and  325  A/m,  respectively  at  /sd 
=  215  pA  and  225  pA.  There  is  current  in  the  negative 
X  direction  in  the  vicinity  of  ymin  at  small  currents  due  to 
the  dominance  of  the  confining  potential.  When 
/sd  =  215  pA  and  225  pA,  however,  E,(y^in)  >  l/tc(ymin)l 
and  no  current  flows  in  the  -x  direction  anywhere  across 
the  sample. 


E 


S' 


Fig.  9.  Current  density  Jt(y)  for  Isd  =  25  p-A,  215  pA,  and  225  pA. 


Fig.  8.  Electric  fields  Edy)  =  -  VVc(y)  and  Et(y)  =  -VVily)  for 
/sd  =  215  pA. 


The  current  I  (y)  for  electrons  moving  in  the  positive 
X  direction  is 


i(y)  = 


My)  dy 


YM) 

Rh 


(27) 


where 


/sd  =  J  7t(y)dy  =  /(y^ax)  + /(ymin) , 

y  min 


(28) 


and 
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We  divide  the  sample  width  into  20  equal  segments  in 
Fig.  10  and  determine  the  percentage  of  current  flowing 
through  each  segment  for  /so  =  25  p.A,  215  p.A,  and  225 
p,A.  We  do  not  show  a  plot  for  Ad  =  0  p-A,  but  A/  would 
be  -0.81  pA  and  +0.81  pA  for  the  left-hand  side  and 
right-hand  side  segments,  respectively,  and  zero  for  the 
other  18  segments  because  /^(rhs)  =  -/c(lhs)  =  0.81  pA. 


Finally,  we  investigate  the  charge-redistribution 
-e8o-(y)  of  the  electrons  in  the  2DEG  in  terms  of  the 
deviation  Scrfy)  in  the  number  of  electrons/cm^  from 
the  average  number  =  5.94X10*Vcm^  on  the  i  =  2 
plateau  at  12.3  T,  where 

ijfi  ^  d^ 


-200  -100  0  100  200 


y  (lini) 


<1 


y  (ym) 


Fig.  10.  Percentage  of  current  flowing  through  20  equal  segments 
across  the  sample  width  for  /sd  =  25  |xA,  215  |xA,  and  225  |xA. 


from  Eq.  (A-5).  Eigure  11  is  a  logarithmic  plot  of 
1 8cr(y )  I  versus  y  for  Ad  =  215  pA.  There  is  an  excess  of 
electrons  on  the  +y  side  of  the  sample,  and  a  depletion 
on  the  -y  side. 

An  assumption  made  in  deriving  Eq.  (A-5)  was  that 
the  charge  density  varies  slowly  across  the  sample,  i.e., 
that  p(y)  =  p(y+dy),  or  8(T(y)  «  n^.  This  assumption 
is  valid  here  because  the  largest  value  of  8cr(y)  occurs 
at  jmax ,  and  is  2  %,  2  %,  and  6  %  of  Ws  at  Ad  =  25  pA, 
215  pA,  and  225  pA,  respectively. 

One  of  the  consequences  of  our  approximate  form  of 
the  charge-redistribution  potential  is  that  the  net  charge 
does  not  vanish  when  the  charge-redistribution  -e  8cr(y ) 
is  integrated  across  the  sample  width.  The  area  under 
the  curves  in  Eig.  11  is  4  %  larger  for  the  +y  side  than 
for  the  -y  side.  Therefore,  there  is  an  unaccounted 
excess  of  electrons;  so  this  is  not  quite  the  actual  shape 
of  the  charge-redistribution  function.  However,  it  is  the 
potential  and  current  distributions  that  are  of  primary 
importance  to  the  transport  properties — not  the  charge- 
redistribution.  The  charge  was  certainly  conserved  in 
the  experiment  of  Eontein  et  al.  [8],  and  yet  their  mea¬ 
sured  potential  distributions  are  symmetrical.  This  fact 
demonstrates  that  the  slight  charge  asymmetry  does  not 
significantly  affect  the  potential  and  current  distribu¬ 
tions. 


The  current  distributions  in  Eig.  10  are  virtually  iden¬ 
tical  between  25  pA  and  215  pA,  even  though  large 
numbers  of  electrons  are  being  excited  into  higher 
Landau  levels  at  215  pA.  The  left  and  right  side  distri¬ 
butions  are  nearly  symmetric.  There  is,  however,  a  sig¬ 
nificant  transfer  of  current  from  the  left-hand  side  seg¬ 
ment  to  the  right-hand  side  segment  at  225  pA.  We  saw 
in  Sec.  5.2  that  no  current  flows  within  60,  55,  and  13 
magnetic  lengths  of  the  sample  side  for  Ad  =  25  pA, 
215  pA,  and  225  pA,  respectively.  Also,  68  %,  70  %, 
and  51  %  of  the  current  is  in  the  19  segments  to  the  left 
of  the  right-hand  side  segment  where  the  edge  channel 
current  would  flow  for  these  three  applied  currents.  The 
current  density  was  negative  in  the  left-hand  side  of 
Eig.  9  at  Ad  =  25  pA  because  electrons  were  flowing  in 
the  -X  direction  at  y^in  -199.554  pm,  but  that  contribu¬ 
tion  to  Al  in  the  left-hand  side  segment  of  Eig.  10  is  so 
small  that  the  net  current  is  positive. 


y  (nm) 


Fig.  11.  Logarithmic  plot  of  the  charge-redistribution  I8a'(jv)l  across 
the  sample  for  /so  =  215  fxA,  where  8(r()7)  is  the  deviation  from  the 
average  number  density  The  ”  region  represents  an  excess  of 
electrons,  the  “-I-”  region  a  depletion  of  electrons. 
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We  could  conserve  the  charge  by  adjusting  the  origin 
slightly  to  the  right  until  the  area  under  the  curves  are 
equal  for  ±y  in  Fig.  11,  and  then  self-consistently  recal¬ 
culating  the  potentials  with  the  new  coordinates.  This 
would  greatly  complicate  the  calculations  however,  and 
with  all  the  approximations  that  have  been  made  in  this 
paper,  and  with  the  excellent  agreement  with  experiment 
[8],  it  seems  unnecessary.  It  may  be  a  consequence  of 
this  charge  nonconservation  problem  that  the  value  of 
ymin  is  inside  the  confining  potential  for  the  case  in 
Table  1  when  /so  =  225  p.A. 

5.4  Sample-Width  Dependence  of  the  Critical 
Current 

Balaban,  Meirav,  and  Shtrikman  [10]  have  found  that 
the  critical  current  for  breakdown  of  the  quantum  Hall 
effect,  la,  scales  logarithmically  with  the  sample  width 
w  for  all  Landau  levels.  We  verify  this  dependence  by: 

(a)  using  the  result  in  Sec.  3.3  that  /c(rhs)  =  -/c(lhs)  = 
0.8 1  p-A  for  the  i  =  2  plateau  at  12.3  T  when  Ad  =  0  p-A; 

(b)  calculating  the  value  of  y^ax  from  Eqs.  (2)  and  (5)  for 
each  value  of  w  at  Ad  =  0  p.A;  (c)  calculating  the  value 
of  G  from  Eq.  (17)  for  each  value  of  w;  (d)  defining  Ar 
as  the  applied  current  Ad  sufficient  to  excite  the  lowest, 
M  =  1,  QUILLS  transitions  [11,25-28];  (e)  assuming 
the  value  of  Ac  (y max)  is  the  same  for  all  values  of  Ad  that 
excite  M  =  I  QUILLS  transitions  (The  value  used  is 
Ac(ymax)  =  5.96X10^  V/m,  obtained  from  Table  1  at 
Ad  =  215  p.A,  w  =  400  p-m,  and£t(ymax)  =  1.1X10'’ V/m. 
This  is  equivalent  to  fixing  the  value  of  A(rhs)  to  be  2.30 
p,A  for  each  value  of  Ar);  (f)  calculating  the  value  of  y^ax 
from  Eq.  (2)  for  each  value  of  w ;  (g)  adjusting  the  value 
of  A  so  that  At(ymax)  =  1.1X10'’  V/m  in  Eqs.  (2),  (17), 
(20),  and  (23);  and  (h)  adjusting  the  value  of  y^in  to  give 
the  correct  Hall  voltage  for  each  current  by  using  Eqs. 
(2),  (15),  (17),  (21),  and  (22). 

The  results  of  Ar  versus  w  are  plotted  in  Eig.  12.  The 
shape  of  the  curve  is  identical  to  the  experimental  data 
of  Balaban  et  al.  [10].  The  scaling  is  very  different, 
however  because  their  critical  currents  are  about  two 
orders  of  magnitude  smaller  than  ours.  We  note  that  the 
experiment  of  Haug,  von  Klitzing,  and  Plog  [29]  tends 
to  agree  with  the  experimental  curve  shapes  of  Balaban 
et  al.  [10],  but  the  experiment  of  Kawaji,  Hirakawa,  and 
Nagata  [30]  found  a  linear,  rather  than  a  logoarithmic, 
dependence  of  Ar  with  w .  Perhaps  this  difference  is  due 
to  nonuniformities  in  the  values  of  the  charge-depletion 
width  A  along  the  sides  of  the  samples,  e.g.,  we  have 
observed  different  values  of  Ar  along  the  lengths  of  some 
of  our  samples.  If  we  assume  that  the  value  of  Vm 
remains  constant  along  a  sample  edge,  allow  A  to  vary 
by  changing  the  average  ionized  donor  density  Ad,  and 
assume  the  ratio  Ac (y max) /At (y max)  remains  constant,  then 


we  find  that  the  critical  current  required  to  excite  M  =  I 
QUILLS  transitions  with  At (y max)  =  LI  X  10'’  V/m 
decreases  when  A  decreases  and  Ad  increases,  i.e.,  the 
steeper  the  confining  potential,  the  smaller  the  critical 
current. 


Fig.  12.  The  critical  current  Ict  versus  sample  width  w.  Refer  to  Sec. 
5.4  for  details  of  the  calculation. 


6.  Conclusions 

We  have  calculated  potential  and  current  distributions 
across  the  width  of  a  GaAs/AlGaAs  heterostructure 
sample  for  applied  currents  between  0  p.A  and  225  p-A, 
using:  (a)  a  quadratic  confining  potential  yc(y)  arising 
from  charge-depletion  regions  along  the  sides  of  the 
sample;  (b)  parameters  for  that  potential  obtained  from 
a  localization  experiment  [9];  (c)  a  logarithmic  charge- 
redistribution  potential  VAy)  of  the  2DEG;  and  (d)  a 
maximum  electric  field  At(ymax)  calculated  from  break¬ 
down  measurements  and  a  QUILLS  model  [11].  Our 
predictions  are  in  excellent  agreement  with  experiments 
[8,10]. 

Referring  to  Table  1,  the  confining  potential  compo¬ 
nent  Ac(ymax)  of  the  electric  field  at  y^ax  contributes 
88  %,  54  %,  and  57  %  to  At(ymax)  at  25  pA,  215  pA,  and 
225  pA,  respectively.  The  maximum  current  density 
■/t(ymax)  is  34  A/m,  85  A/m,  and  325  A/m,  respectively  at 
these  three  currents.  A  significant  amount  of  current  is 
distributed  within  the  sample  interior.  Eor  example,  A  is 
99  %,  99  %,  and  83  %  of  Ad,  respectively  at  these  three 
currents.  We  predict  the  current  to  be  much  farther  from 
the  sides  of  the  sample  than  in  other  models,  e.g.,  no 
current  flows  within  60,  55,  and  13  magnetic  lengths  of 
the  sample  side  for  these  currents.  It  would  require  a 
lateral  resolution  of  about  0. 1  pm,  rather  than  the  70  pm 
resolution  of  Eontein  et  al.  [8],  to  verify  this  result. 
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7.  Appendix  A.  Derivation  of  Eqs.  (10) 
and  (11) 

We  first  derive  Eq.  (10)  of  Sec.  4.  The  Lorentz  force 
exerted  on  the  conducting  electrons  causes  deviations 
-e8cr(y')  from  the  average  surface  charge  density 
-ews  =  -  ie^B/h  of  the  2DEG  at  each  point  y'  across  the 
sample  width.  Consider  a  strip  of  this  redistributed 
charge  of  width  dy located  in  the  x-y  plane  at  position 
y’  and  pointing  in  the  x  direction,  with  a  charge/length 
A{y')  =  -e8(T(y')dy'.  Gauss’s  law  is  then  used  to  obtain 
the  electric  field  dE  (y )  at  some  point  y  in  the  x-y  plane, 

due  to  one  of  these  line  charges:  ksq^E  ■  dS  =  q,  and 

thus  KeodE(y)27rly'-ylL;,  =  A{y')L^  for  a  cylindrical 
Gaussian  surface  of  radius  \y'-y\  and  length  L^.  The 
potential  dy(y)  of  this  line  charge  is 


dV(y)  =  -  d£(y)dy  =  - 


A{y’)  2  r  1 

liTKeoW  J  2ly-yl/w 

y-w/2 


dy 


=  Aizlin^ly^ 
2ttkso  w 


(A-1) 


Fig.  A-1.  Geometry  factors  for  the  calculations  in  Appendix  A. 


If  the  volume  charge  density  is  slowly  varying,  then 
p(y)  =  P(y  +  where  ph  =  -ecr  and  8ph  =  -ehcr. 
Thus,  from  Eqs.  (A-3)  and  (A-4), 


or 


8t7(y)dy  =  ^  [£(y)-£(y  +  dy)], 

(x>ctS 


8o-(y)  =  - 


o-ly)  d£’(y) 
w^B  dy 


<r{y)  d^V(y) 
(d^B  dy^ 


The  total  potential,  when  summed  over  all  the  line 
charges,  is  Eq.  (10) 

wH 

=  /§-(/)  In 

-w/2 

(A-2) 

Now  we  derive  Eq.  (11)  of  Sec.  4.  We  found  in  Eq.  (3) 
of  Ref.  [11]  that  the  center  of  mass  coordinate  yo  of  each 
state  undergoing  cycloidal  motion  in  a  Landau  level  is 


-  i/-yi 

W 


dy. 


But  cr(y)  =  ieB/h,  so  we  obtain  Eq.  (11) 

j.  ,  ,  _  d^V(y)  _  im*  dV(y) 
“  hwe  dy"  “  hB  dy" 


(A-5) 


Equation  (11)  considers  the  charge-redistribution 
-e8cr(y)  due  to  the  second  derivative  of  the  charge- 
redistribution  potential  V^iy).  The  charge-redistribution 
that  we  calculated  in  Sec.  5.3  depends  also  on  the 
second  derivative  of  those  regions  of  the  confining 
potential  yc(y )  which  differ  from  the  y^ax  and  ymm  values 
at  /sD  =  0  p-A. 


yo  = 


E(y)  . 

m,B 


fik^  . 

eB  ’ 


(A-3) 


so  the  states  move  to  the  right  as  E(y')  increases.  The 
total  charge  8Q  transferred  into  the  volume  hhdy,  out¬ 
lined  with  solid  lines  in  Eig.  A-1,  is 

8Q  =  8phhdy  =  kh[Ayo(y) p  (y) 


-Ayo(y  H-  dy)p (y  H-  dy)].  (A-4) 


8.  Appendix  B.  Eigenvalue  Equation 

We  saw  in  Eq.  (12)  of  Sec.  4  that 


K(y)  =  -  f  I  ^  [K(y')]  In 


-  l/-yl 

W 


dy',  (B-1) 


where  the  characteristic  length  ^  is 


2TTKS(,h(Dc 


(B-2) 
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MacDonald,  Rice,  and  Brinkman  [19]  have  pointed  out 
that  Eq.  (B-1)  has  the  form  of  an  eigenvalue  equation 
®[^(y)]  =  CF(y),  where 


0(y)  =  - 


w/2 


[2  ,  ,  ,1 

J  5?=' 

-  ly  -yi 

Lw  J 

-w/2 


is  an  operator  acting  on  the  function  F(y)  =  V,  (y'),  and 
C  =  1/^  is  a  constant. 

The  function  F(y)  =  V,{y')  should  have  a  second 
derivative  that  satisfies  Eq.  (B-1).  Our  function,  from 
Eq.  (15)  of  Sec.  4.1,  is 


K(y)  =  -^ 


In 


v-Sm 


In 


y’  +  w/2 
y'  -  w/2 


(B-3) 


where  the  cut-off  is  =  ^  -  Jmax  , 


dy 


-,K(y')  =  - 


4Rh 

2 


In 


-  W 

[(w/2)^-(/)^]’ 


and 


d" 

dy'^ 


Vr(y') 


FRn 

2 


In 


2wy’ 

[{wl2f~{y'ff 


(B-4) 


We  can  see  if  the  potential  given  by  Eq.  (B-3)  is  a 
valid  solution  to  the  eigenvalue  expression  Eq.  (B-1)  by 
substituting  Eq.  (B-4)  into  Eq.  (B-1)  and  integrating 
only  between  the  limits  -  S^in  to  S^ax  because  8cr(y')  is 
zero  beyond  these  two  cut-off  values.  Surprisingly,  we 
obtain  nearly  exact  solutions  to  the  eigenvalue  ex¬ 
pression  at  all  values  of  y  when  ^is  less  than  2.0  pm  and 
5max  =  Smin  =  This  choice  of  S^ax  =  S^m  =  ^  for  the 
cut-off  values  was  used  by  Beenakker  and  van  Houten 
[22].  Our  value  of  ^  is  3.3  nm,  so  we  would  be  well 
within  this  exact  range  if  5n,ax  =  S^m  =  The  values  of 
Smax  do  not  equal  5n,i„  however,  and  are  much  larger  than 
the  value  of  Also,  we  have  a  confining  potential, 
ydy),  parts  of  which  should  be  included  in  Eq.  (B-1). 
It  would  be  interesting  to  see  how  well  the  values  of 
Vt(y)  obtained  in  Eq.  (B-4)  agree  with  the  values 
obtained  in  the  eigenvalue  expression  Eq.  (B-1)  when 
using  the  values  of  quantities  obtained  in  Table  1. 
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